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c5 : ABSTRACT 

0^ I The bulk of the quiet solar photosphere is thought to be significantly magne- 

I tized, due to the ubiquitous presence of a tangled magnetic field at subresolution 

' scales with an average strength (i?)~100 G. This conclusion was reached through 

C/^ ■ detailed three-dimensional (3D) radiative transfer modeling of the Hanle effect in 

^ ■ the Sr I 4607 A line, using the microturbulent field approximation and assuming 
that the shape of the probability density function of the magnetic field strength is 
exponential. Here we relax both approximations by modeling the observed scat- 

Iri • tering polarization in terms of the Hanle effect produced by the magnetic field of 

I— I' a 3D photospheric model resulting from a (state-of-the-art) magneto-convection 

^ ■ simulation with surface dynamo action. We show that the scattering polarization 

^ ■ amplitudes observed in the Sr I 4607 A line can be explained only after enhanc- 

ly^ ' ing the magnetic strength of the photospheric model by a sizable scaling factor, 

^ ■ F^IO, which implies (S)~130 G in the upper photosphere. We argue also that 

^ ■ in order to explain both the Hanle depolarization of the Sr I 4607 A line and the 

O ' Zeeman signals observed in Fe I lines we need to introduce a height-dependent 

^ ■ scaling factor, such that the ensuing {B)^160 G in the low photosphere and 

j> ■ (-B)~130 G in the upper photosphere. 
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1. Introduction 

The importance of the Hanle effect as a diagnostic tool of the quiet Sun magnetism 
cannot be overestimated. This magnetic-field-induced modification of the linear polarization 
due to scattering in spectral lines is even sensitive to the presence of magnetic fields that 
are randomly oriented at scales too small to be resolved (Stenflo 1982). Applying novel 
diagnostic techniques based on the Hanle effect, Trujillo Bueno et al. (2004) reached the 
following conclusions (see also Trujillo Bueno et al. 2006): 

• The quiet solar photosphere is permeated by a small-scale tangled magnetic field, whose 
average strength at heights h^^SOO km above the visible solar surface is (S)f=:il30 G 
when no distinction is made between granular and intergranular regions. 

• The downward-moving intergranular lane plasma is pervaded by relatively strong tan- 
gled magnetic fields at subresolution scales, with (B) > 200 G. This conclusion suggests 
that most of the fiux and magnetic energy reside on still unresolved scales. 

• In the upper photosphere the ensuing energy flux estimated using the typical value of 
lkms~^ for the convective velocity (thinking in rising magnetic loops) or the Alfven 
speed (thinking in Alfven waves generated by magnetic reconnection) turns out to be 
substantially larger than that required to balance the radiative energy losses from the 
quiet chromosphere. 

The first conclusion was obtained by contrasting low resolution observations of the 
scattering polarization signals of the Sr I 4607 A line with three-dimensional (3D) radiative 
transfer calculations in a realistic 3D model of the quiet solar photosphere obtained from the 
hydrodynamical (HD) simulations of Asplund et al. (2000) (hereafter the HD model). The 
second conclusion resulted from combining the information provided by the strontium line 
with that we obtained through the application of the Hanle effect line ratio technique for 
C2 lines explained carefully in Trujillo Bueno et al. (2006). It is however important to note 
that the above-mentioned HD model is unmagnetized and that, therefore, such conclusions 
are based on the following hypotheses for the quiet Sun magnetic field that produces Hanle 
depolarization: (1) the magnetic field is tangled at scales smaller than the mean free path 
of the line-center photons, with an isotropic distribution of directions and (2) the shape of 
the probability density function (PDF), describing the fraction of the quiet Sun photosphere 
occupied by magnetic fields of strength B, is exponential (PDF(i?) = -^e"^/^^^). 

The main aim of this paper is to determine the strength of the magnetization of the 
quiet solar photosphere without making use of such two approximations. To this end, we 
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investigate the Hanle effect of the Sr I 4607 A hne in a 3D photospheric model resulting 
from the magneto- hydrodynamic (MHD) simulations of Vogler & Schiissler (2007), which 
show a complex small-scale magnetic field that results from dynamo amplification of a weak 
seed field. These authors went a step further than Cattaneo (1999) by demonstrating that 
a realistic fiow topology of stratified, compressible, surface convection without enforced re- 
circulation is capable of local dynamo action near the solar surface. It is interesting to note 
that the small-scale "hidden magnetic field" inferred via the Hanle effect lends support to 
the notion of a turbulent solar surface dynamo playing a significant role in the magnetism 
of the quiet Sun (e.g., Pietarila Graham et al. 2009; see also the reviews by Lites 2009, and 
by Sanchez Almeida & Martinez Gonzalez 2011). 

2. Formulation of the Hanle-efFect problem in 3D magneto-convection models 

We calculated the number density of Sr I atoms at each grid point of a 3D magneto- 
hydrodynamic model of the solar photosphere resulting from a non-gray radiative transfer 
version of the Run-C simulation of Vogler & Schiissler (2007) (hereafter, the MHD model). 
To this end, we solved the standard 3D non-LTE problem without polarization including all 
the allowed radiative and coUisional transitions between the 15 bound levels of a realistic 
model of Sr I and the ensuing ionizing transitions to the ground level of Sr II. With the 
resulting number density of Sr I atoms fixed, we solved the full 3D radiative transfer problem 
of the scattering polarization in the Sr I 4607 A line by applying the complete frequency 
redistribution theory of resonance line polarization (e.g., Landi DeglTnnocenti & Landolfi 
2004), which is suitable for modeling the fractional linear polarization observed in weak 
resonance lines like that of Sr I at 4607A (see Sampoorna et al. 2010; and references therein). 

The Sr I line at 4607 A results from a single ^Sq — ^Pi resonance line transition with 
Ji — and Ju — 1 (the total angular momentum of its lower and upper levels, respec- 
tively). Although the ensuing Einstein coefficient for spontaneous emission is very large (i.e., 
A„;Ri2xl0^s~^) the line originates in the solar photosphere because the abundance of stron- 
tium is relatively low (see Figure 1). Scattering polarization in the Sr I 4607 A hne is caused 
only by the atomic polarization (population imbalances and quantum coherences) of its up- 
per level, which results from anisotropic radiation pumping in the solar photosphere. The 
lower level cannot be polarized and the relevant transfer equations for modeling the available 
low resolution observations of the fractional linear polarization signals arc -^X = X — Sx, 
with X the Stokes parameter I, Q or U at the frequency and direction of propagation under 
consideration {dr — —rjids defines the monochromatic optical path along the ray, rij being 
the absorption coefficient and s the geometrical distance). 
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The source function components are Si — r S}^^ + {1 — r)Bi,, Sq — r Sq^^ and Su — 
r 5"^°^, where B^, is the Planck function and r = k^^i// {ni(l>u + f<-c) (with ki the hne-integrated 
opacity, Kc the continuum opacity and (pi^ the normahzed absorption profile). The expres- 
sions of the line source function components arc identical to those used by Trujillo Bueno 
& Shchukina (2007), which are functions of the inclination and azimuth of the ray under 
consideration and of the following position-dependent quantities (with Q = for K = and 
Q = 0, 1,2 for X = 2): 5J = ^ J^j^+i ^Q (^^^^ Pq multipolar components of the 
upper- level density matrix normalized to the overall population of the ground level). Note 
that Sq (with Q — 1,2) are complex quantities; thus, we will use Sq and Sq to denote, re- 
spectively, the real and imaginary parts of Sq. The values of such quantities at each spatial 
grid point have to be found by solving the following coupled set of equations (Manso Sainz 
& Trujillo Bueno 1999): 



S^ = (l-e)J° + eB,, 
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where e = Cui/(A„; -|- Cui) (with C^; the rate of inelastic coUisions with electrons, in s"-*^), 
Fu = 8.79 X lO^QuB/Aui (with Qu the Lande factor of the upper level and B the magnetic 
strength in gauss), 5^'^^ = D^'^^Aui (with D^'^'^ the upper-level depolarizing rate due to elastic 
collisions with neutral hydrogen atoms), and w^^q = 1. In these expressions, Jq (with Q = 
for K = and Q = 0, 1, 2 for A' = 2) arc the spherical components of the radiation field tensor 
(Landi Degl'Innocenti & Landolfi 2004), which are integrals over the frequency and direction 
of the Stokes parameters weighted by the normalized absorption profile. Thus, Jq quantifies 
the familiar mean intensity of the incident radiation, Jq its anisotropy, while Jq and Jq 
(with Q — 1,2) denote, respectively, the real and imaginary parts of Jq (which measure the 
breaking of the axial symmetry of the incident radiation field) . The Mij-coefficients depend 
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on the inclination {9b) of the local magnetic field vector with respect to the vertical Z-axis 
and on its azimuth (xb)- 

The depolarizing rate, 5^^^ = D^'^yA^i, is proportional to the neutral hydrogen number 
density (see Eq. 33 in Faurobert-SchoU et al. 1995). Such formula for calculating 5^^) 
results from an accurate quantum-mechanical derivation, and we point out that it leads 
to values similar to those we have obtained through the semi-classical theory of Anstee & 
O'Mara (1995). The excellent agreement between the 5^^^ values obtained from the quantum 
and semi-classical approaches is a clear indication of the reliability of the 5*^^^ values of our 
modeling, because the Sr I 4607 A line is a triplet type transition with spin S=0 (i.e., it is a 
spectral line for which the classical description holds). 

Of particular interest is the fact that the magnetic strength needed to reach any given 
depolarization is the larger the greater 6^"^^ This can be seen in the expression of the 
critical magnetic field strength for the onset of the Hanle effect, which can be defined as the 
strength (Be) for which the Hanle depolarization reduces the zero-field scattering polarization 
amplitude by about a factor of two. Eq. (A16) of Trujillo Bueno & Manso Sainz (1999) 
implies that 

B,^{1 + 5^'^)Bh, (3) 

where Bh = 1.137x10^^ Aui/qu is the critical Hanle field in the collisionless regime (i.e., 
Bh^23 G for the Sr I 4607 A line). We point out that in the atmospheric region of formation 
of the Sr I 4607 A line (i.e., approximately between 250 km and 350 km for line of sights with 
0.6 > > 0.1 in the MHD model) the horizontally averaged Be value varies between 60 G and 
40 G, approximately, with (Be) decreasing with height because {S^"^^) decreases with height 
in the solar atmosphere (following the density decrease) . More or less similar B,. values are 
found in one-dimensional semi-empirical models. It is also important to note that a lower 
limit for the Hanle saturation field of the Sr I 4607 A line is ^satur — 200 G ~ 10 Bh, since 
full saturation is reached only for B>300 G. 

3. Results 

The results reported here were obtained through the self-consistent numerical solution 
of the equations, applying the methods mentioned in Trujillo Bueno & Shchukina (2007). 
In a first step we solved the scattering line polarization problem of the Sr I 4607 A line in 
the MHD surface dynamo model, but imposing 5 = G at each spatial grid point. The 
resulting Q/I and U/I line-center signals are similar to those shown for three line-of-sights 
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in figure 1 of Trujillo Bueno & Slicliukina (2007), wliicli were obtained by solving tlie same 
3D radiative transfer problem in tlie above-mentioned HD model. For the moment, the 
scattering polarization observations of the Sr I 4607 A line that have been published lack 
spatial resolution (see Figure 2). We point out that when the /, Q and U profiles we have 
calculated in the MHD model are spatially averaged over scales significantly larger than that 
of the solar granulation pattern we obtain U/I^O and the Q/I line-center amplitudes given 
by the black open circles of Figure 2. Remarkably, these reference scattering polarization 
amplitudes are similar to those computed by Trujillo Bueno et al. (2004) in the HD model 
(see the green open circles) @ The small but noticeable differences between the two curves 
are due to the fact that the thermodynamic structure of the MHD model is not as accurate 
as in the HD model, with Jq/Jq and 1/(1 + 5*^^^) being slightly larger in the MHD model 
(in the photospheric region where the Q/I profile of the strontium line originates, whose 
amplitude is approximately proportional to the product of such two quantities). 

The black dashed line shows the Q/I amplitudes that we obtain when taking into 
account the Hanle depolarization produced by the actual magnetic field of the MHD model, 
whose magnetic Reynolds number is i?in~2600. Not surprisingly, the Hanle depolarization 
produced by the model's magnetic field is too small to explain the observations. This is 
because the mean field strength at a height of about 300 km in the MHD model is only 
(B) ?al5 G (i.e., an order of magnitude smaller than the (B) value inferred by Trujillo 
Bueno et al. 2004). Clearly, the level of small-scale magnetic activity of the MHD surface 
dynamo model is significantly weaker than that of the real quiet Sun photosphere. 

In order to estimate the mean field strength of the quiet Sun photosphere we scaled 
the magnetic strength at each spatial grid point by a constant factor, F, till fitting the 
observations of Figure 2. We have done this keeping the magnetic field topology of the MHD 
model fixed. We find F^12, which implies (B) ^150 G at a height of 300 km above the 
"visible surface" in the MHD model (i.e., an average field strength similar to that inferred 
by Trujillo Bueno et al. 2004). Note that this result has now been reached without using 
the two hypotheses mentioned in §1. 

Curiously enough, a similar scaling factor is also found when one follows Pietarila Gra- 
ham et al. (2009) and assumes that the ratio between the magnetic strengths of the real solar 
photospheric plasma and of the original MHD surface dynamo model is equal to the square 
root of the ratio of the corresponding magnetic Reynolds numbers, with i?m(Sun) f»3xl0^ 
and i?m (model) ~ 2600. As shown by the dotted line of Figure 3, to scale the field strength of 



^For this reason we find again (B) wl30 G when the same hypotheses mention in §1 are used to infer {B) 
from calculations in the B = G version of the MHD model. 
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the MHD model with a constant scaling factor F^12 implies an unrealistic height- variation 
of (B) . The ensuing large magnetic field strengths in the region of formation of the Zeeman 
signals of the Fe I lines at about 15650 A and 6302 A (e.g., around 60 km in the MHD model 
for the Fe I 6302.5 A line; see Figure 1) would produce synthetic Stokes profiles in contra- 
diction with those observed. In fact, F = 12 is significantly larger than the tentative scaling 
factor 3 needed by Danilovic et al. (2010) for explaining the histograms of the polarization 
signals produced by the Zeeman effect in the 6301.5 A and 6302.5 A Fe I lines. On the 
other hand, F = 3 is too low a value for explaining the scattering polarization observations 
of Figure 2 (see the dashed-dotted line). Therefore, in order to be able to fit the observa- 
tional data of Figure 2 without destroying the possibility of explaining the Stokes profiles 
observed in the Fe I lines, we have introduced a height-dependent scaling factor which leads 
to the solid line of Figure 3 for the height variation of the mean field strength. As shown by 
the black solid line of Figure 2, the magnetic field that results from our height-dependent 
scaling factor produces a Hanle depolarization in the Sr I 4607 A line that explains well 
the observation^. It is also important to mention that a significant fraction of the model's 
granular plasma that contributes to the scattering polarization of the C2 lines mentioned in 
§1 is magnetized with {B) ~10 G. 



4. Conclusions 

We have solved the radiative transfer problem of resonance polarization and the Hanle 
effect of the Sr I 4607 A line in a 3D model of the quiet solar photosphere resulting from the 
magneto-convection simulations with surface dynamo action of Vogler & Schiissler (2007). 
We find that the level of magnetic activity in their surface dynamo model is too low for 
explaining the scattering polarization observations of the Sr I 4607 A line. The observed 
linear polarization amplitudes can however be explained after multiplying each grid-point 
magnetic strength by a scaling factor F^12. 

A scaling factor F^12 is about three or four times larger than that required to explain 
the polarization induced by the Zeeman effect in the Fe I lines around 6302 A. This discrep- 
ancy between the scaling factors needed to explain the Hanle signals of the Sr I 4607 A line 
and the Zeeman signals of the Fe I lines disappears if we take into account that the Zeeman 
signals originate a few hundred kilometers deeper than the Hanle signals (see Figure 1). In 



^The small discrepancy between the black solid line and the data points of Figure 2 can be explained 
by the fact that the thermal and density structure of the MHD model is not as accurate as that of the HD 
model. 
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other words, with a height-dependent scahng factor that imphes (B) fa 160 G in the low pho- 
tosphere and (B) 130 G in the upper photosphere it is possible to explain approximately 
both types of observations. 

Our investigation of the Hanle effect in the above-mentioned surface dynamo model 
reinforces the conclusions of Trujillo Bueno et al. (2004) summarized in §1. In particular, 
the horizontal variation of the scaled magnetic field needed to fit the observations indicates 
that most of the flux and magnetic energy reside on still unresolved scales in the intergran- 
ular plasma. Moreover, the fact that the magnetic energy density carried by this "hidden" 
magnetic flcid is a signiflcant fraction of the local kinetic energy density and that the scat- 
tering polarization amplitudes of the Sr I line and of the C2 lines do not seem to be clearly 
modulated by the solar cycle (Trujillo Bueno et al. 2004; Kleint et al. 2010), support 
the suggestion that a small-scale turbulent surface dynamo plays a significant role in the 
magnetism of the "quiet" Sun. 
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Fig. 1. — The regions of formation, in the MHD model, of the (Stokes Q) Hanle signals 
of the Sr I 4607 A line and of the (Stokes V) Zeeman signals of the Fe I 6302.5 A line. 
The colored background provides information on the vertical and horizontal variation of 
the model's magnetic strength, distinguishing between positive and negative values of the 
field's vertical component. The dotted line around /i~0 km indicates the height where the 
continuum optical depth at A4607 is unity. The solid lines indicate the heights where the 
line-center optical depth of the Sr I 4607 A line is unity along line of sights with /i = 0.6 
(thick line) and /i = 0.1 (thin line). The dashed-dotted line shows the heights where the 
optical depth of the Fe I 6302.5 A line is unity along line of sights with // = 1, calculated at 
the wavelengths where the modulus of the Stokes V profiles have their maximum value. 
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Fig. 2. — Center-to-limb variation of the (spatially averaged) Q/I scattering amplitudes 
of the photospheric line of Sr I at 4607 A. The data correspond to various observations 
taken by several authors during a minimum and a maximum of the solar activity cycle (see 
references in Trujillo Bueno et al. 2004). The two green lines show scattering polarization 
amplitudes calculated in the HD model, without including any magnetic field (green open 
circles) and including the Hanle depolarization of a microturbulent field with an exponential 
PDF characterized by a mean field strength (i?) = 130G (green dashed line). The black 
lines show scattering polarization amplitudes calculated in the MHD model, neglecting its 
magnetic field (black open circles) and taking into account the Hanle depolarization of the 
model's magnetic field with scaling factors F = 1 (black dashed line) and F = 3 (black 
dashed-dotted line). As shown by the black solid line, the observations can be approximately 
fitted by multiplying each grid-point magnetic strength by a height-dependent factor F{h) 
(which implies the height- variation of (B) given by the solid line of Figure 3). Practically 
the same black solid line is obtained with a constant scaling factor F = 12 (which implies the 
unrealistic height- variation of (B) given by the dotted line of Figure 3). The red solid line 
shows the calculated scattering polarization amplitudes when imposing B > i?satur = 300 G 
at each grid point in the MHD model. The fact that this red line coincides with 1/5 of the 
zero-field Q/I amplitudes indicates that the microturbulent field approximation is indeed a 
suitable one. 
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Fig. 3. — Various height-variations of the mean field strength corresponding to four different 
scahng factors F of the magnetic strength of the MHD modeL Dashed-hne: F = 1 (i.e., 
as in the original MHD model). Dashed-dotted line: F = 3 (i.e., as suggested by the 
Zeeman polarization investigation of Danilovic et al. 2010). Dotted line: F = 12. Solid 
fine: a height-dependent scaling factor that explains the scattering polarization observations 
of the Sr I 4607 A fine without destroying the possibility of explaining simultaneously the 
polarization induced by the Zeeman effect in the Fe I fines. The arrows at /iRi60 km and 
/i~300 km indicate the approximate atmospheric heights around which the observed Zeeman 
and Hanlc signals are produced, respectively. Note from the solid and dashed lines that 
F{h = 60 km) w 4 while F{h = 300 km) w 10. 



